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Abstract In order to better understand the solar genesis of interplanetary 
magnetic clouds (MCs) we model the magnetic and topological properties of 
four large eruptive solar flares and relate them to observations. We use the three- 



dimensional Minimum Current Corona model (Longcope, 1996 ) and observations 
of pre-flare photospheric magnetic field and flare ribbons to derive values of 
reconnected magnetic flux, flare energy, flux rope helicity and orientation of 
the flux rope poloidal field. We compare model predictions of those quantities to 
flare and MC observations and within the estimated uncertainties of the methods 
used find the following. The predicted model reconnection fluxes are equal to or 
lower than the reconnection fluxes inferred from the observed ribbon motions. 
Both observed and model reconnection fluxes match the MC poloidal fluxes. The 
predicted flux rope helicities match the MC helicities. The predicted free energies 
lie between the observed energies and the estimated total flare luminosities. The 
direction of the leading edge of the MC's poloidal field is aligned with the poloidal 
field of the flux rope in the AR rather than the global dipole field. These findings 
compel us to believe that magnetic clouds associated with these four solar flares 
are formed by low-corona magnetic reconnection during the eruption, rather 
than eruption of pre-existing structures in the corona or formation in the upper 
corona with participation of the global magnetic field. We also note that since all 
four flares occurred in active regions without significant pre-flare flux emergence 
and cancellation, the energy and helicity we find are stored by shearing and 
rotating motions, which are sufficient to account for the observed radiative flare 
energy and MC helicity. 

Keywords: Flares, Relation to Magnetic Field; Helicity, Magnetic; Flares, Mod- 
els 



1. Introduction 

Coronal mass ejections (CMEs) expel plasma, magnetic flux and helicity from 
the Sun into the interplanetary medium. At 1 AU in the interplanetary medium 
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CMEs appear as interplanetary coronal mass ejections (ICMEs). At least one 
third (Gosling, 1990) or perhaps a larger fraction (Webb et al., 2000) of the 
ICMEs observed in situ are magnetic clouds (MCs) ( |Burlaga et al., 19"81~ ), 
coherent "flux-rope" structures characterized by low proton temperature and 
strong magnetic field whose direction typically rotates smoothly as they pass 
the satellite. 

MCs originate from eruptions of both quiescent filaments and active regions 
(ARs). The 3D magnetic models and geomagnetic relationships are better un- 



derstood for filament eruptions than for ARs ( 


Marubashi, 1986 


Bothmer and 


Schwenn, 1998 


Zhao and Hoeksema, 1998 


Yurchyshyn et al, 2001 


). However 


the most geoeffective MCs originate from ARs ( 


Gopalswamy et al, 2010 


). In 



this paper we focus exclusively on the latter. 

Comparison of the properties of MCs with those of their related ARs clarifies 
our understanding of both domains. Assuming MCs to be twisted flux ropes in 
magnetic equilibrium, several authors have succeeded in inferring global proper- 
ties such as MC axis orientation, net magnetic flux, and magnetic helicity (see 
review by Demoulin (2008)). For a sample of twelve MCs Leamon et al. (2004) 
found that the percentage of MC poloidal flux relative to unsigned vertical AR 
flux varied widely, from 1% to 300%. For one MC |Luoni et al\ ( |2005[ ) did a 
similar study and found a factor 10 times lower flux in the MC than the AR, in 
agreement with other previous studies (Demoulin et al., 2002 Green et al., 2002). 
More recently, for a sample of nine MCs, Qiu et al. ( 2007 1 found that the MC 
poloidal flux matches the "observed" reconnection flux, i.e. reconnection flux in 
the two- ribbon flare associated with it, and the toroidal flux is a fraction of the 
reconnection flux. The jQiu et al. (20071 results may be interpreted as evidence 
of formation of the helical structure of magnetic flux ropes by reconnection, in 
the course of which magnetic flux, as well as helicity, is transported into the flux 
rope. 

The other quantity that is very useful for relating MCs to their associated 
flares is magnetic helicity, which describes how sheared and twisted the magnetic 



field is compared to its lowest energy state (Berger, 1999 Demoulin and Pariat, 



2007). Since helicity is approximately conserved in the solar atmosphere and the 



heliosphere (Berger and Field, 1984), it is a very powerful quantity for linking 
solar and interplanetary phenomena. For six ARs Nindos, Zhang, and Zhang] 
(2003) found that photosphcric helicity injection in the whole AR is comparable 
with the MC helicity. However, it is worth remembering that this approach is 
simplified, since the liftoff of the flux rope does not remove all of the helicity 
available in the AR (Mackay and van Ballegooijen, 2006 Gibson and Fan, 2008 ). 



Mandrini et al. ( 2005 1 and Luoni et al. ( 2005 ) compared, respectively, the helicity 



released from a very small AR and a very large AR, with the helicity of their 
associated MC. They found a very good agreement in the values (small AR with 
small MC, and large AR with large MC), despite a difference of 3 orders of 
magnitude between the smaller and the larger events. 

There exist two basic ideas about the solar origin of magnetic clouds: MCs 
are formed either globally or locally. In the global picture, the MC topology is 
defined by the overall dipolar magnetic field of the Sun (Crooker, 2000). In this 
case, the field lines of the helmet streamer belt become the outermost coils of the 
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MC through reconnection behind the CME as it lifts off. Hence the leading field 
direction of magnetic cloud tends to follow that of the large-scale solar dipole, 



reversing at solar maximum (Mulligan, Russell, and Luhmann, 1998 Li et al., 



2010 Bothmcr and Schwenn, 1998). In the local picture, on the other hand, the 



flux rope is formed within the AR and its properties are defined by properties of 
the AR. We can categorize the "local" models into two sub-classess. In the first, 
the magnetic flux rope emerges from beneath the photosphere into the corona 



(Low, 1994 Chen, 1989 Fan and Gibson, 2004 Leka et al., 1996 Abbett and 



Fisher, 2003). In this scenario the flux ropes formed may maintain stability 



for a relatively long time prior to the explosive loss of equilibrium (Forbes 



and Priest, 1995 Lin, Raymond, and van Ballegooijen, 2004) or a breakout 



type reconnection that opens up the over lying flux rope of opposite polarities 
(Antiochos, Devore, and Klimchuk, 19991. Such a flux rope is therefore pre- 
existing before its expulsion into interplanetary space. In the second case the 
flux rope is formed in situ by magnetic reconnection. The magnetic reconnection 
suddenly re-organizes the field configuration in favor of expulsion of the "in situ" 
formed magnetic flux rope out of the solar atmosphere. The results of |Qiu et al.\ 
( 2007 ) support this case. Qiu et al. ( 2007 ) found that the reconnection flux from 
observations of flare ribbon evolution is greater than toroidal flux of the MC but 
comparable and proportional to its poloidal flux, regardless of the presence of 
filament eruption. Their conclusion agrees with the inference from the study by 



Leamon et al. (2004), although through a very different approach. 



Our working hypothesis is that MCs associated with the ARs originate from 
the ejection of locally in-situ formed flux ropes. In this case shearing and rota- 
tion of the photosphere magnetic flux concentration before the flare lead to the 
buildup of magnetic stress which is removed during the flare by reconnection. 
As a result a magnetic flux rope is formed and erupts, producing a MC. 

To test our hypothesis we apply a quantitative non-potential self-consistent 



model, the Minimum Current Corona (MCC) model (Longcope, 1996 Longcope, 



2001 ), to predict the properties of the in-situ formed flux rope in four two-ribbon 



flares. Using the MCC model with SOHO/MDI magnetogram sequences we 
construct a three-dimensional model of the pre-flare magnetic field topology and 
make quantitative predictions of the amount of magnetic flux that reconnects in 
the flare, the magnetic self-helicity of the flux rope created, and the minimum en- 
ergy release the topological change would yield. We then compare the predicted 
flare helicity and energy to MC helicity and flare energy, inferred from fitting 
the magnetic cloud (Wind/ ACE) and GOES observations correspondingly. We 
compare the predicted reconnected flux to the amount of photospheric flux swept 
up by the flare ribbons using TRACE 1600 A data and the poloidal MC flux 
inferred from fitting the magnetic cloud observations. We find that for the four 
studied flares our results support, from the point of view of flux, energy and 
helicity, the scenario in which the MC progenitor is a helical flux rope formed in 
situ by magnetic reconnection in the low corona immediately before its expulsion 
into interplanetary space. We also find that MC topology is defined by the local 
AR structure rather than the overall dipolar magnetic field of the Sun in the 
events studied. 

This paper is organized as follows: in Section [2] we describe the methods and 
uncertainties of our analysis. In Section[3]we describe the four flares studied, the 
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Table 1. Flare and AR Properties (see jj3|. Number is the NOAA number of the AR 
associated with the flare, &ar, in units of 10 22 Mx, is the AR's unsigned magnetic 
flux and Har, in units of 10 42 Mx 2 , is helicity injected into the AR during the 
magnetogram sequence starting at to and ending at tji are . 



i 




Flare 






Active Region 


M-gram sequence 




Date 


Time 


Class 


Number 


®AR 


H A r 




1 


05/13 2005 


16:57 


M8 


10759 


2.0 


-12 ± 1.2 


05/11 23:59 
















05/13 16:03 


2 


11/07 2004 


16:06 


X2 


10696 


2.1 


-15 ± 1.5 


11/06 00:03 
















11/07 16:03 


3 


07/14 2000 


10:03 


X6 


09077 


3.4 


-27 ± 2.7 


07/12 14:27 
















07/14 09:36 


4 


10/28 2003 


11:10 


X17 


10486 


7.5 


-140 ± 14. 


10/26 12:00 
















10/28 10:00 



flux and helicity of the ARs in which the flares occurred, and the magnetogram 
sequence during the buildup time. In Section [4] we discuss our results, and in 
Section [5] summarize our conclusions. 



2. Methods: Calculating Reconnection Flux, Energy and Helicity 



In this section we describe the methods that we use to ({ 2.1 1 predict the recon- 



nection flux, energy and helicity from SOHO/MDI magnetogram sequences and 



the Minimum Current Corona model and ([2.2 1 determine observed values of 



these quantities from GOES, TRACE, ACE, and WIND observations. 



2.1. Minimum Current Corona Model 



The key improvement of our study relative to Leamon et al. (2004) and Qiu 



et al. (2007) is the use of the Minimum Current Corona model, which allows us 



to estimate the energy and helicity associated with the in-situ formed flux rope 



(Longcope, 1996 


Longcope, 2001) 


of the CSHKP ( 


Carmichael, 1964 


Pneuman, 1976) 


two-ribbon flare 



Sturrock, 19681 Hirayama, 1974 Kopp and 



shearing of an AR along its polarity inversion line (PIL) to build up stress. After 
this pre-flare stress buildup, the MCC model quantifies the result of eliminating 
some or all of the stress and creating a twisted flux rope overlying the AR, 
through magnetic reconnection. 

To describe the evolution of the pre-flare photospheric motions that lead to 



stress build-up we use a sequence of SOI/MDI full-disk magnetograms (Scherrer 



et al., 1995). As the starting point we take t , right after the end of a large 
flare, which we call the zero-flare. We assume that at to the AR's magnetic field 
becomes fully relaxed. As the ending time we take tfi are , right before the time 
when the flare of study occurred but avoiding artifacts associated with the onset 
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Table 2. Flare physical properties: MCC model predictions vs. observations: pre- 
dicted (& r ,MCc) and inferred from observations (& r ,ribbon) reconnection fluxes and 
MC poloidal fluxes ($ p ,mc), predicted (Smcc) and observed (Egoes) energy values, 
predicted (Hmcc) and observed (Hmc) helicity values (see 



ribbon 



p,MC 



10 21 Mx 10 21 Mx 10 21 Mx 

1 2.8 ±0.4 4.1 ±0.4 6.3 ±4.2 

2 5.4 ±0.8 4.8 ± 0.5 5.25 ± 3.5 

3 6.0 ±0.9 12.8 ±3 9.9 ± 6.6 

4 15.0 ± 2.6 23 ±7 18.0 ± 12.0 



£-mcc Egoes H M qc H M c 

10 31 ergs 10 31 ergs 10 42 Mx 2 10 42 Mx 2 

1.0±0.3 -7.0±1.2 

3.1 ±0.6 -7.5 ±5.0 

6.4±1.8 -5.0±0.6 

2.0 ±0.1 -8.3 ±5.5 

9.1 ±2.6 -20.1 ±3.6 

10.1 ±2.1 -22.5 ±15.0 

18.0 ±5.2 -48.0 ±8.6 

13.6 ±0.6 -45.0 ±30.0 



of the flare brightening ( Qiu and Gary, 2003 ) . To achieve the maximum energy 
release, the field reconnecting during the flare would need to relax to its potential 
state, hence we assume the field to be potential at tfi are . As a result, we form 
a sequence of magnetograms, which covers At hours of stress build-up prior to 
the flare (see Table [l}. 

For quantitative analysis of the pre-flare magnetic field we divide each mag- 
netogram into a set of unipolar partitions and then into unipolar magnetic 
charges (e.g. see partitioned magnetogram in Appendix, Figure|4|. Firstly, for all 
successive pairs of magnetograms we derive a local correlation tracking (LCT) 
velocity field ( November and Simon, 1988 Chae, 2001 ) and then group pixels into 
individual partitions that have persistent labels. In the second step we represent 
each magnetic partition with a magnetic point charge (or magnetic point source) 
which has the flux of the partition and is located at its center of flux. Finally, 
using the LCT velocity field we calculate the helicity injected by the motions of 
the magnetic point charges of the whole AR, (Har, see Tableland Longcope 



et al. (2007)). We make sure that the amount of helicity injected by the motions 



of the continuous photospheric partitions matches the helicity injected by the 
motions of the magnetic point charges. Their equality gives us confidence that 
the centroid motions of the point charges accurately capture helicity injection. 
Computing the vector potential via the Fourier approach of Chae (2001), as we 
choose to do, results in the higher values (10%) in the helicity flux compared to 
that from the approach by Pariat, Demoulin, and Berger ( 2005 1 ( Chae, 2007 ). In 
addition, the LCT method that we use yields systematically lower values than 
the DAVE velocity inversion algorithm with a difference in helicity flux of less 
than ~ 10% (Welsch et al., 20071. Those two effects result in an uncertainty of 
10% in the Har value which we take into account. 

The MCC model characterizes the changes in the pre-flare magnetic field 
purely in terms of the changes in the magnetic domains, volumes of field lines 
connecting pairs of opposite point charges. Replacing each partition with a single 
magnetic point charge as we chose to do results in values of domain fluxes that 
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are only slightly different from the actual domain fluxes (Longcope, Barnes and 



Beveridge, 2009). As magnetic charges move, the magnetic field, first relaxed 
by the zero-flare, becomes increasingly stressed and hence non-potential. Un- 
der the assumption that no reconnection, flux emergence or cancellation occur 
between the zero-flare and the flare of interest, the domain fluxes could not 
have changed. (Note, that for our analysis we selected only flares associated 
with ARs with no significant flux emergence or cancellation during the period 
between to and tfi are ). To provide both the domain flux conservation and the 
increasing field non-potentiality, the MCC model includes currents only on the 
intersections between the domain boundaries, called separators. In this way the 
lack of reconnection leads to storage of free magnetic energy, energy above that 
of the potential field, which could then be released by reconnection in the flare. 
To achieve the maximum energy release, the field inside the domains associated 
with the flare (flaring domains) would need to relax to its potential state. Thus 
to find the reconnection flux we first need to find the flaring domains and then 
calculate the changes in their domain flux from t to tfi are . More specifically, 
we first overlay the magnetic point charges rotated to the time of the TRACE 
1600 A flare ribbons onto the ribbons image and then use a Monte Carlo method 
(Barnes, Longcope, and Leka, 2005) to find the fluxes of the flaring domains at 
to and tfi are . Finally, we separately sum up the absolute values of all the positive 
and negative changes in the domain fluxes to calculate the model reconnection 
flux ($r.MCC, see Table |2]). This is the model estimate of the net flux transfer 
that must occur in the two-ribbon flare through the flare reconnection. 

To find the flaring separators we find the topology of the magnetic field at 
tfiare and select those separators that connect nulls that are located on the flare 
ribbons. Through the MCC model, the changes in the domain fluxes under those 
flaring separators allow us to calculate current, free energy and helicity liberated 
on each separator (for a detailed description of the method see Longcope ( 1996) 



and Appendix B of Kazachenko et al. (2009)). The total model energy (Smcc) 



released during the flare is a sum of energies released at each flaring separator. 
It is a lower bound on the energy stored by the pre-flare motions, since MCC 
model uses the point charge representation and hence applies a smaller number 
of constraints than point-for-point line-tying. It can be shown that the energy 
stored by ideal, line-tied, quasi-static evolution will always exceed the energy of 
the corresponding flux constrained equilibria ( Longcope and Magara, 2004 1 . The 
total mutual helicity injected on all flaring separators is a sum of the helicities 
injected on each flaring separator. However the liftoff of the flux rope does not 



remove all of the helicity available in the flux rope ( Mackay and van Ballegooijen, 



2006). For an MHD-simulated eruption Gibson and Fan (2008) found that 41% 



of the helicity is lost with the escaping rope, while 59% remains. For simplicity, 
we assume that 50% of the total mutual helicity from the MCC model ends 
up as self helicity of the flux rope created by reconnection: Hmcc = 
Finally wc note that the MCC model depends on the way we partition the 
magnetogram sequence (Beveridge and Longcope, 2006). We experimented with 
different values of saddle points in the partitioning and apodizing windows in the 
LCT and found that this contributes an uncertainty in the MCC reconnection 
flux, MCC energy and MCC helicity that we include in Table [2j 
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2.2. Flare And MC Observations 



For comparison to the predictions of the MCC model, we must infer values 
of reconnection flux $ r ,ribbon from the observations, MC poloidal flux $ p ,MCi 
energy £ go es and helicity Hmc (see Table [2]) . 

To infer values of reconnection flux $ r ribbon from the observations we use flare 



ribbon motion (Poletto and Kopp, 1986 Fletcher and Hudson, 2001) observed 
m 1600 A images from TRACE. To find the total magnetic flux swept out by a 
moving ribbon, we count all pixels that brightened during any period of the flare 
and then integrate the unsigned magnetic flux encompassed by the entire area 
taking into account the height of the ribbon's formation, a « 20% correction 
(Qiu et al, 2007). The uncertainties in the $ r ,ribbon are estimated by artificial 



misalignment between the MDI and TRACE data, ribbon edge uncertainty and 
inclusion of transient non-ribbon features with the ribbon areas. To quantify 
the misalignment contribution we perform a set of trials whereby magnetogram 
and 1600 A images are offset by up to 2 MDI pixels. To find the uncertainty 
due to ribbon edge identification we perform the calculation for different ribbon- 
edge cutoff values ranging from 6 to 10 times the background intensity. We also 
compare the MCC reconnection flux & r ,Aicc to the poloidal MC flux $ Pi mc 
derived from fits to the in-situ MC ACE/ Wind observations using the Grad- 
Shafranov reconstruction method (Hu and Sonnerup, 2001 ). As Qiu et al. (2007) 
showed from observations, <& r ,ribbon ~ &p,MC- Hence if the MCC model captures 
the reconnection flux correctly, & r ,MCC should match $ p ,mc unless reconnection 
of the ICME with the ambient solar wind makes an important contribution 



(Dasso et al, 2006). 



During the flare the magnetic free energy that has been slowly stored by 
photospheric motions £mcc is rapidly released by reconnection and then dissi- 
pated. We estimate energy losses not only due to radiation (£ r 



as 



Kazachenko 

et al.\ Kazachenko et al. ( 2009| 2010) did, but also due to conductive cooling 
(£ c ) and the enthalpy flux (£ e nt), which in some numerical cases is as large as 
radiative energy losses (Bradshaw and Cargill, 2010 ). Since it is not clear whether 



the source for the CME kinetic energy is the magnetic free energy stored in the 
active region and not the energy stored e.g. in the interplanetary current sheet, 
we neglect the energy carried away by the CME. 

To quantify the three components of £qo es = £ r + £c + £ent (see Table [3| we 
use GOES analysis software in SolarSoft and the observed GOES X-ray fluxes in 
the two channels (1-8 A and 0.5-4 A). Those provide an estimate of the plasma 
temperature T and emission measure EM — n^V, where n e is the electron 
density and V is the emitting volume. Radiative energy losses £ r depend on the 
emission measure, temperature and composition of emitting plasma. We find 
their magnitude using the temperature dependent Mewe radiative loss function 
(Mewe, Gronenschild, and van den Oord, 1985). To calculate the conductive 



energy losses £ c we integrate the conductive energy loss rate to the chromosphere 
Pcond = Uth/Tcond where Uth is thermal energy content of the plasma Uth = 
2>n e kTV = 3kTy/EM x V and T con d is the cooling time scale 



7~cond 



3fcn e (L/2) 2 
k T 5 /2 



(1) 
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Table 3. Observed energy budget (in 10 ergs): radiative losses (£ r ), conductive 
losses (£ c ), enthalpy fluxes (Sent), total energy (£goes = £ c + £ r + £ent and 
estimated value for flare luminosity £fl ~ (3.15 ± 1.05) x £goes- The p redic ted 
model energy (£mcc) is given for comparison with the observations (see f 2.2 and 



i L (Mm) S r 



ent 



Zgoes 



£fl 



Smcc 



145 ±31 1.0 

43 ± 20 0.9 

151 ± 50 2.5 

107 ±18 5.3 



0.45 ±0.05 1.6 ±0.6 

0.4 ±0.1 0.7 ±0.1 

2.6 ±0.2 4.9 ±1.9 

2.75 ± 1.35 5.5 ±0.7 



3.1 ±0.6 10.3 ±5.1 

2.0 ±0.1 6.4 ±2.4 

10.1 ±2.1 34.0 ±17.2 

13.6 ±0.6 43.5 ±16.1 



1.0 ±0.3 
6.4 ± 1.8 

9.1 ± 2.6 
18.0 ±5.2 



for a loop of full length L, Boltzmann constant k and Spitzer conductivity kq ~ 
10 -6 (Longcope et al., 2010). We quantify the volume of the emitting material V 



by assuming that £ c « £ r at late times, as they should be in a static equilibrium 



(Rosner, Tucker, and Vaiana, 1978 Vesecky, Antiochos, and Underwood, 1979). 



From the volume V and emission measure EM we derive the electron density 
n e = WEM/V. For the loop length L we use the distribution of the lengths 
of the flaring separators (with the energy weights) whose geometrical properties 
are found from the coronal magnetic topology at tfi are (|2.1|. The mean and 



standard deviation of the lengths of the flaring separators yield the mean and 
the standard deviation of the values of £ c . Finally, we estimate the enthalpy 
flux £ ent using model calculations by Bradshaw and Cargill (2010). From Tables 



1 and 2 in the paper by Bradshaw and Cargill (2010) and the loop lengths of 



the flaring separators, we first derive a coefficient S which describes the ratio 
between the radiative cooling and the enthalpy flux time scales and then the 
enthalpy flux itself. 

We get an additional idea for the value of the uncertainty in £qo es by compar- 
ing it to the flare luminosity (FL, £fl) from the Total Irradiance Monitor (TIM) 
on the Solar Radiation and Climate Experiment (SORCE). Unfortunately, FL 
measurements are not available for the four flares studied here. However FLs 
have been measured for four other large (>X10) solar flares (see Table 2 in 



Woods, Kopp, and Chamber lin (2006)), for which we may calculate £goes- For 



these four flares we find that the FLs are approximately two to four times larger 
than the £ G oes- We use this scaling range to limit our energy estimates from 
above (see Table [3}: £ F l ~ (3.15 ± 1.05) x £ GO es- 

Finally, we compare the model MCC flux rope helicity with the helicity of 
the magnetic cloud associated with the flare, Hmc- We calculate Hmc applying 



the Grad-Shafranov method ( |Hu and Sonnerup, 2001] ) to the ACE/Wind MC 
observations. There are several uncertainties and limitations in the determina- 
tion of Hmc, which are as well applicable to MC poloidal flux calculations. 
First, the inferred value of Hmc is model-dependent: e.g. within the cylindrical 
hypothesis, force-free and non-force-free models give helicities values that differ 



by up to 30% (Dasso et al., 2003 Dasso et al., 20061. However, this variation 



remains small compared to the variation of helicity values computed for different 
MCs (Gulisano et al., 2005). Second, the MC boundaries can be defined by 
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several criteria, which do not always agree. This introduces an uncertainty on 
the magnetic helicity which can be comparable to the uncertainty obtained with 
different models (Dasso et ai, 2006). Finally, the distribution of the twist along 
the flux rope, as well as the length of the flux rope are generally not known. So far 
only in one case the length of the flux rope Lmc — 2.5AU has been determined 



precisely from impulsive electron events and solar type III radio bursts (Larson 



et ai, 1997). We take the value of 0.5 AU as the lower limit of Lmc (DeVore, 



2000) and 2.5 AU as the upper limit of Lmc (Larson et al., 1997). Such choices 



of the lower and upper limits of Lmc would change poloidal MC flux and MC 
helicity to vary between roughly half and twice the measured values. 



3. Data: Flares Studied 

We apply the methods described in the previous section to four large eruptive 
flares (Table [l]). This number of events is limited by several necessary flare 
selection criteria. Firstly, we selected only events which have good observations 
of both the flare and the MC. Secondly, except for the May 13 2005 flare, we 
selected only ARs where two successive flares larger than M-class were present, in 
order to make plausible our assumption of initial relaxation of the AR's magnetic 
field to potential state. Thirdly, both the flare of study and the zero-flare should 
happen no farther than 40° from the central meridian so that the stress build- 
up could be observed. Finally, we selected only flares associated with ARs with 
no significant flux emergence or cancellation during the period between to and 

t flare- 

Our topological analysis using the MCC model has been executed previously 
for three of four flares: M8 flare on May 13 2005 (Kazachenko 2009), X2 flare 



on November 7 2004 (Longcope et al., 2007), X17 flare on October 28 2003 



(Kazachenko et al., 2010). The results of the MCC analysis for the X5.7 flare 
on July 14 2000 are described in this paper for the first time (see Appendix). In 
Table[l]we list the flare number (i) in this work, date, time and X-ray class of each 
flare; the NOAA number of the AR (AR) associated with the flare, AR's unsigned 
magnetic flux ($>ar) and helicity injected into the AR during the magnetogram 
sequence (Har); start (to) and end time (tfi are ) of the magnetogram sequence. 
The flares are sorted by X-ray class. In Table[2]we compare MCC model predicted 
physical properties with the observations: predicted and inferred from the ribbon 
motions reconnection fluxes and MC poloidal fluxes, predicted and observed 
from the GOES observations energy values, predicted and observed from the 
Wind/ ACE observations helicities. Finally in Table [3] we detail the observed 
energy budget for each flare. 

The first flare l isted in Table [j] is the M8 flare that occurred on May 13 2005 



in NOAA 10759 (Kazachenko et ai, 2009 Yurchyshyn et ai, 2006 Liu et ai, 



2007 Jing et al., 2007 Liu, Zhang, and Zhang, 2008). NOAA 10759 has a large 



positive sunspot which contains more than a half of the total positive flux of the 
AR and rotates with the rate of 0.85° ± 0.13° per hour during 40 hr before the 



flare (Kazachenko et al., 2009). Such fast rotation along with the fact that the 



spin helicity flux is proportional to the magnetic flux squared makes the effect 
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of sunspot rotation dominant in the helicity budget of the whole AR. As for the 
flare itself, the rotation of the sunspot produced three times more energy and 
magnetic helicity than in the hypothetical case in which the sunspot does not 
rotate; the inclusion of sunspot rotation in the analysis brings the model into 
substantially better agreement with GOES and interplanetary magnetic cloud 
observations. Rotation is energetically important in the flare and alone can store 
sufficient energy to power this M8 flare. 

The second flare in Table [l] is the X2 flare on November 7 2004 ( Longcope 



et al., 2007). The start time was plausibly taken to be that of an M9.3 flare 



which occurred 40 hr before the flare of interest. The MCC model predicts a 
value of the flux needed to be reconnected in the flare that compares favorably 
with the flux swept up by the flare ribbons. The MCC model places a lower 
bound on the energy stored by the 40-hour buildup shearing motions that is at 
least three times larger than the observed energy losses. The helicity assigned to 
the flux rope that is assumed in the model is comparable to the magnetic cloud 
helicity. Note that our estimate for Hmcc m Table [2] is higher than the one 
in Longcope et al. (2007) (see Table II in Longcope et aL|(2007)): we estimate 



Hmcc as a sum of the helicities over all eight flaring separators, while |Longcope| 



et al. (2007) took a sum over only the three most energetic separators. 



The third flare in Tabic [T] is the X6 flare on July 14 2000 ( |Lepping et~aL, 



2001] |Yurchyshyn et al, 2001] |Fletcher and Hudson, 200T] |Masuda, Kosugi, and 



Hudson, 2001 Sol. Physics special edition, 20011. Our analysis of this flare is 



described in this paper for the first time (see Appendix). We take as the zero- 
flare an XI. 9 flare around 48 hr before the flare of interest. We use the MCC 
model to find that the released energy is comparable to the observed energy 
losses. The amount of flux reconnected during the flare according to the model 
is at least one and a half times smaller than the reconnection flux observed with 
TRACE. The model estimate for the helicity is comparable with the helicity 
from the MC observations. No sunspot rotation is associated with the pre-flare 
evolution. 

The fourth event in Table [l] the X17 Halloween flare (Yurchyshyn, Hu, and 



|Abramenko, 2005[ |Regnier and Priest, 200?] |Schrijver et al., 2006[ |Mandrini| 



et al, 2006 Lynch et al, 2005 Zhang, Liu, and Zhang, 2008), occurred in an 



AR with a fast-rotating sunspot. We find that the MCC reconnection flux is 
consistent with the reconnection flux inferred from the observations. We find 
that the sunspot rotation increases the total AR helicity by 50%. However 
in contrast to the flare on May 13 2005, where rotation is dominant in the 
energetics, rotation increases the free energy and flux rope helicity of this flare 
by only w 10%. Shearing motions alone store sufficient energy and helicity to 
account for the flare energetics and ICME helicity content within their observa- 
tional uncertainties. Thus this flare demonstrates that the relative importance 
of shearing and rotation in this flare depends critically on their location within 
the parent AR topology. 
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4. Results: MCC Model Predictions Versus Observations 

The main global property that describes the flare's reconnection is the amount of 
magnetic flux that participates. Figure [T] shows predicted ($ r ,MCc) and observed 
(^ r ., ribbon) reconnection fluxes for each event, the AR average unsigned magnetic 
flux ($>ar) and the poloidal MC flux (& p ,mc)- The first thing we notice is 
that the fraction of the AR magnetic flux that is observed to reconnect during 
the four flares ranges from 18% to 49%. Secondly, in the second and fourth 
flares the predicted reconnection flux matches the reconnection flux inferred 
from the observations, while in the first and third flares the highest probable 
value of the MCC reconnection flux is lower than the lowest probable value 
of the observed reconnection flux by 13% and 29% correspondingly. The lower 
model reconnection flux is likely due to additional reconnections not accounted 
for in the model. That means that the MCC model captured a lower limit of 
the amount of magnetic flux that has reconnected in these flares and hence 
the lower limit on the amount of energy released. Finally, in all four cases the 
value of poloidal MC flux matches both the observed and model reconnection 
fluxes, although the uncertainties in <& p ^mc due to the unknown MC length are 
quite large. According to the CSHKP model, on which the MCC model builds, 
reconnection contributes solely to the incremental poloidal component of the 
flux-rope flux. Therefore, the derived agreement between the poloidal MC flux 
and the reconnection fluxes means that the flux rope is formed in situ. 

The MCC model gives a lower limit of the free magnetic energy released 
in each flare (Longcope, 2001). In Figure [2] we compare the predicted MCC 



model free energy {Emcc, diamonds) with the observed time-integrated sum of 
radiative and conductive energy losses and the enthalpy flux (Sgoes, stars); we 
also show the estimated flare luminosity (£fl, blue squares). Figure [2] indicates 
that for the third and fourth flares the predicted energy Emcc matches the 
observed energy Sgoes, while for the second flare the predicted energy is around 
three times larger than Sgoes- In all four cases the flare luminosity is higher 
than both £ go es and Emcc ■ This is not surprising, since the MCC model uses 
a point charge representation rather than line tying and yields a lower limit 
on the free energy released in the flare. Summarizing, within the uncertainties, 
for three flares i = 2,3,4 the predicted free energy lies between the observed 
estimate of released energy and the estimated flare luminosity. Only for the May 
13 2005 flare is the model energy lower than the observed estimate. Note, that 
this flare is the only flare that did not have a zero-flare at to. Since rotation is 
the dominant source of helicity injection in this flare and the rotation rate was 
around zero before to, we believe that our analysis plausibly captures the major 
source of helicity injection in this flare. Nevertheless, the absence of a zero-flare 
indicates that there might have been additional energy storage before to- Hence, 
our estimate is a lower limit to the reconnection flux and the magnetic energy 
of this M8 flare. 

One must keep in mind that one of the basic assumptions of the MCC model 
is the potentiality of the magnetic field after the zero-flare. |Su, Golub, and Van] 



Ballegooijen (20071 analyzed TRACE observations of 50 X- and M-class two- 



ribbon flares and found that 86% of these flares show a general decrease in 
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Figure 1. Predicted and observed magnetic flux values for the four events. The horizontal 
axis shows the maximum X-ray flux of each flare. For details see Tables [T] and [2] For discussion 



see 







the shear angle between the main polarity-inversion line and pairs of conjugate 
bright ribbon kernels. They interpreted this as a relaxation of the field towards 
a more potential state because of the eruption that carries helicity/current with 
it, but one can readily argue that a similar decrease in shear angle would be 
seen if sequentially higher, less-sheared post-flare loops light up with time as the 
loops cool after reconnection. These results are consequently ambiguous: they 
may show a decrease in shear, or they may reflect that flares generally do not 
release all available energy and part of the flux-rope configuration remains. In 
other words, the MCC model potentiality assumption may mean that additional 
energy and reconnection flux is stored before the zero-flare. 

Because magnetic helicity is approximately conserved in the corona, even 
in the presence of reconnection, it is instructive to compare Hmcc to Hmc- 
Figure [3] shows the relationship between the predicted (HmcCi diamonds) and 
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Figure 2. Predicted and observed energy values for the four events. For details see Table [3] 
For discussion see jj4] 



observed (Hmc , stars) values of helicity. The blue +s show the amount of the 
helicity for the whole AR (Har)- In all four cases we get MCC model helicities 
that are of comparable magnitude and same sign as the observed MC values 
and are smaller than the helicity of the whole AR. Our analysis shows that 
prcflarc motions contribute enough stress to account for observed helicity values, 
however more accurate estimate of the MC length is required in order to lower 
the error bars in the MC helicities and improve our understanding of MC/flare 
relationship. It is interesting that the most energetic of the four flares (i = 4, 
which happened in the southern hemisphere) had the helicity sign opposite to the 
hemispheric helicity preference ( Pevtsov and Balasubramaniam, 2003k . Hence its 



sign cannot be predicted from the global solar properties, but only from a case 
study like this. 

The properties of the magnetic field in a MC are determined by the initial 
conditions of the eruption which we derive with the MCC model as well as by how 
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Figure 3. Predicted and observed helicity values for the four events. For details see Table [T] 
and Table [2] For discussion see j|4] 



the MC interacts with the interplanetary medium during its travel toward the 
Earth. Above we found a consistency between the MC flux and predicted model 
and observed flare reconnection fluxes, the MC and predicted flare flux rope 
helicities, observed and predicted flare energy releases. The agreement between 
those supports our local in-situ formed flux rope hypothesis. 

One more quantity that is frequently compared between the MC and AR 



flux rope is the direction of the poloidal field. Li et al. (2010 1 found that the 
poloidal field of MCs with low axis inclination relative to the ecliptic (~ 40% of 
all MCs) has a solar cycle dependence. They note that during the solar minima, 
the orientation of the leading edge of the MC is predictable: it is the same as 
the solar dipole field. However during the maximum and the declining phases, 
when most of the geoffective MCs happen, both (north and south) orientations 
are present, although the global dipole field orientation of the beginning of the 
cycle dominates. 
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It is instructive to consider how the four flares of our study relate to these 
results. Assuming that the poloidal field in the ejected flux rope is oriented 
along the flaring separators, we define its orientation relative to the ecliptic 
plane using the North-South classification: north (30° < 9 < 90°, B z > 0, in the 
solar ecliptic coordinate system) or south (—90° < 6 < —30°, B z < 0). We then 
determine the orientation of the leading edge of the MC poloidal field using the 
same North-South classification and compare two quantities: the orientation of 
the poloidal field in the active region and the orientation of the leading edge 
of the MC poloidal field. We find that the MC produced by the Bastille day 
flare during the solar maximum has a south oriented leading MC poloidal field, 
same as both the remnant weak dipole orientation and poloidal field orientation 
of the flux rope at the sun. In contrast, the flares which occurred during the 
declining phase, on May 13 2005 and November 7 2004, produced magnetic 
clouds with south oriented leading MC poloidal fields, opposite to the direction 
of the global dipole field, but same as the poloidal field orientation predicted 
for a flux rope in the modeled AR. Finally the MC produced by the Halloween 
flare laid perpendicular to the ecliptic plane and thus was not relevant to the 
observed Li et al. (2010) rule; however a good agreement was also found between 
the directions of the poloidal field in the MC and in the source AR ( |Yurchyshyn,| 



Hu, and Abramenko, 2005). 



Summarizing the above, although there is a tendency for ARs to follow the 
dipole field orientation during the solar minimum (Li et al., 2010), during the 
solar maximum and the declining phase, when the largest MCs occur, the local 
AR field is important. We find that for the four studied large events the direction 
of the leading MC poloidal field is consistent with the poloidal field orientation in 
the AR rather than to the global dipole field in agreement with |Leamon, C anficld, 



and Pevtsov (2002). This implies that the magnetic clouds associated with large 



ARs inherit the properties of the AR rather than those of the global dipole field, 
as a result of reconnection in the active region rather than with the surrounding 
dipole field. Although here we compare the poloidal post-flare arcade field with 
the poloidal MC field, this supports the conclusion by Yurchyshyn et al. ( |2007 ), 
who found that 64% of CMEs are oriented within 45° to the MC axes (MC 
toroidal field) and 70% of CMEs are oriented within 10° to the toroidal field of 
EUV post-flare arcades (Yurchyshyn, Abramenko, and Tripathi, 2009). In other 
words, despite the fact that CME flux ropes may interact significantly with the 
ambient solar wind (Dasso et al., 2006) or other flux ropes (IGopalswamy et al., 



2001), a significant group of MCs reflects the magnetic field orientation of the 



source regions in the low corona. 



5. Conclusions 

The main purpose of this study is to understand the mechanism of the CME 
flux rope formation and its relationship with the MC. Notably, we use the 
Minimum Current Corona model (Longcope, 19961 which, using the pre-flare 
motions of photospheric magnetic fields and flare ribbon observations, quantifies 
the reconnection flux, energy and helicity budget of the flare. We apply this 
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model to four major eruptive solar flares that produced MCs and compare the 
predicted flux rope properties to the observations. 

We compare model predictions to observations of four quantities: the pre- 
dicted model reconnection fluxes to the MC poloidal fluxes and ribbon motion 
reconnection fluxes, the predicted flux rope helicities to the MC helicities, the 
predicted released energies to the total radiative/conductive energy losses plus 
the enthalpy fluxes, the direction of the magnetic field in the AR arcade to the 
direction of the leading edge of MC poloidal field. 

Our comparison reveals the following. The predicted reconnection fluxes matchj 
the reconnection fluxes inferred from the observations for the November 7 2004 
and Halloween flares. For the May 13 2005 and Bastille day flares the minimum 
probable differences between the predicted and observed reconnection fluxes are 
13% and 29% correspondingly. In all four cases the values of poloidal MC fluxes 
match both the observed and the model reconnection fluxes. The predicted flux 
rope helicities match the MC helicities. For three flares of study the predicted 
free energies lie between the observed energy losses (radiative and conductive 
energy losses plus the enthalpy fluxes) and the flare luminosities. Only for the 
flare on May 13 2005, the predicted free energy is one third of the observed 
estimate. We relate this mismatch to the fact that May 13 2005 flare was the 
only event without a zero-flare, hence additional energy might have been stored 
before to. Finally, we find that in all four cases the direction of the leading MC 
poloidal field is consistent with the poloidal component of the local AR arcade 
field, whereas in two cases the MC poloidal field orientation is opposite to that 
of the global solar dipole. 

These findings compel us to believe that magnetic clouds associated with 
these four eruptive solar flares are formed by low-corona magnetic reconnection 
during the eruption, rather than eruption of preexisting structures in the corona 
or formation in the upper corona by the global field. Our findings support the 
conclusions of Qiu et al. (2007) and Leamon et al. (2004), although through 



a very different approach: while Qiu et al. (2007) and Leamon et al. (2004) 



inferred the solar flux rope properties only from observations, we infer them 
from both the MCC model and the observations. Using the pre-flare magnetic 
field evolution and the MCC model, we find that we are able to predict the 
observed reconnection fluxes within a 29% uncertainty and the observed MC 
poloidal flux and helicity values within the MC length uncertainty. For the flares 
associated with zero-flares we are able to estimate a lower limit for the free 
magnetic energy. We note that, since all four flares occurred in ARs without 
significant pre-flare flux emergence/cancellation, the flux/energy /helicity we find 
is stored by shearing and rotating motions, which is sufficient to account for 
observed energy and MC flux and helicity. 

Our work brings up several interesting questions that require further explo- 
ration. Firstly, the results of this paper are based on only a small number of 
events, which are similar in that all have a large radiative signature. Hence a 
study of observations of a class of events with small radiative signature would 
be challenging: smaller flares that are nevertheless associated with major CMEs 
( Aschwanden,Wuelser, and Nitta, 2009). If the MCC model is valid, it should 
be able to explain both the energy and helicity content of flare/CME events 
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whose flare energy output is disproportionately small. Secondly, in this paper, we 
estimated the final flare energy as a sum of energy losses by radiation, conduction 
and enthalpy, neglecting the energy carried away by the CME. To our knowledge, 
no systematic empirical study of the source of energy for the CME has yet been 
conducted, so it is unknown what proportion of its energy budget arises from 
the AR. From the limited sample, Ravindra and Howard (20101 found that 
a 50% contribution may be a reasonable first-order approximation. Including 
CME energy losses into our flare analysis would help us understand how much 
of the CME energy arises from the active region and may lead to a greater 
understanding of the onset mechanism for CMEs. 



Appendix 

The X5.7 Bastille Day flare occurred on 2000 July 14 2000 at 10:03 UT in 
NOAA 9077. Our magnetic field data describing the evolution of the magnetic 
field before this flare consist of a sequence of SOI/MDI full-disk magnetograms 
(2", level 1.8) starting at t = 2000 July 12 14:27 UT, after the XI. 9 flare 
(2000 July 12 10:18 UT), and ending at t flare =July 14 09:36 UT, 27 minutes 
before the Bastille day X5.7 flare. Thus we form a sequence of 28 low-noise 
magnetograms with a 96-minute cadence, which cover 43 hours of the stress 
buildup prior to the X5.7 flare on 2000 July 14 10:03 UT. Firstly, for all suc- 
cessive pairs of magnetograms we use a Gaussian apodizing window of 7" to 
derive a local correlation tracking (LCT) velocity. We then take a magnetogram 
at tfi are and group pixels, exceeding a threshold B t hr = 45 Gauss downhill 
from each local maximum, into individual partitions. We combine partitions 
by eliminating any boundary whose saddle point is less than 350 Gauss below 
either maximum it separates. Each partition is assigned a unique label which 
maintains through the sequence by using the LCT velocity pattern. Figure [4] 
shows the spatial distribution of these partitions at tfi are . For expediting the 
assessment of the field's connectivity we represent each magnetic partition with 
a magnetic point charge which contains the magnetic flux of the whole partition 
concentrated in the partition's centroid. We find that the magnetic field is well 
balanced at t// ore ($+($_) = 3.3(— 3.5) x 10 22 Mx) and exhibits no significant 
emergence/cancellation during the 43 hours of pre-flare stress buildup time. From 
the LCT velocity and magnetic field in each point we find the flux of relative 
helicity into the corona to be H AR = -(27 ± 2.7) x 10 42 Mx 2 , no significant spin 
helicity content (rotation) has been detected. 

To find the model estimate of the reconnection flux we determine the mag- 
netic point charges associated with the flare using the flare UV observations by 
TRACE 1600 A. Figure|5]shows a superposition of the elements of the topological 
skeleton at tfi are onto the UV flare ribbon image. The spines (red solid lines) that 
are associated with ribbons form the footprint of a combination of separatrices 
which overlay the flaring domains. The overlay suggests that the northern ribbon 
is associated with the spines connecting flaring point sources P08, Pll, P20, 
P01, P02, P16, P04, P03; and the southern ribbon is associated with the spines 
connecting N20, N13, N02, Nil, N04, N07, N06, N09, N10, N12, N03, N01, N18, 
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Figure 4. Positive (P) and negative (N) polarity partitions for NOAA 9077 on July 14 09:36 
UT. The gray-scale magnetogram shows magnetic field scaled from -1000G to 1000G. The 
partitions are outlined and the centroids are denoted by +'s and x's (positive and negative 
respectively). Axes are labeled in arc-seconds from disk center. 



N19, N15, N05. Field lines connecting the pairs of opposite point charges listed 
above form a set of flaring domains. The amount of flux that those domains 
exchanged, the model reconnection flux, is & r ,MCC = (6.0 ±0.9) x f0 21 Mx, fifty 
percent smaller than the lower value of the observed reconnection flux from the 
flare ribbon evolution ($ r ,ribbon = 12.8 ± 3 x 10 21 Mx). 

From the set of flaring point charges and nulls lying between them we find 
twenty six flaring separators (see Figure [6]). The total free energy and helicity 
output on those separators is £mcc = (9.1 ± 2.6) x 10 31 ergs and Hmcc = 
-(20.1 ± 3.6) x 10 42 Mx 2 . However, out of the 26 flaring separators 90% of the 
total free energy is contained in separators originating in nulls B08 and Bll 
which lie between poles P01, P02 and P16 (note the most red separators in 
Figure [6]). Moreover, 64% of the total flare free energy is partitioned between 
three separators: A19/B11, A20/B11, A20/B08. According to the MCC model 
the poles associated with those nulls (P01, P02, P16 and N10, N01, N03, N18) 
indicate the locations of the largest free energy release. Figure [6] indicates that 
the brightest observed loops in TRACE 1600 A are the loops connecting point 
charges N12, N10, N01 and N18 with P02, consistent with the results of the 
MCC model presented above. 
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Figure 5. TRACE 1600 A image, plotted as reverse gray scale, with elements of the topo- 
logical skeleton superimposed. The skeleton calculated for July 14 09:36 UT is projected onto 
the sky after its tangent plane has been rotated to the time of the TRACE observations (10:33 
UT). Positive and negative sources are indicated by +'s and x's respectively. The triangles 
represent the labeled null points. The red curved line segments show spine lines associated 
with the reconnecting domains. Axes are in arc-seconds from disk center. 
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